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SUMMARY 


An iterative method has been developed for designing wing 
section contours corresponding to a prescribed subcritical distri- 
bution of pressure. The calculations are initialized by using a 
surface panel method to analyze a baseline wing or wing-fuselage 
configuration. A first-order expansion to the baseline panel 
method equations is then used to calculate a matrix containing the 
partial derivative of potential at each control point with respect 
to each unknown geometry parameter. In every iteration cycle, the 
matrix is used both to calculate the geometry perturbation and to 
analyze the perturbed geometry. The distribution of potential on 
the perturbed geometry is established by simple linear extrapola- 
tion from the baseline solution. The extrapolated potential is 
converted to velocity by numerical differentiation, and velocity 
is converted to pressure by Bernoulli's equation. Not only is the 
accuracy of the approach good for very large perturbations, but 
the computing cost of each complete iteration cycle is sub- 
stantially less than one analysis solution by a conventional panel 
method. Example design solutions are presented to demonstrate 
that the method is accurate, numerically stable, and converges in 
only three to five iterations. 


1 . INTRODUCTION 


The surface panel approach has proved to be very successful 
for the analysis of subcritical, potential flow around geometri- 
cally complex aircraft configurations (references 1-10). It is 
well established that the better formulated panel methods consis- 
tently predict accurate wing pressure distributions, including the 
effect of strong fuselage-nacelle interference. A typical example 
is illustrated in Figure 1. Several investigators who have been 
aware of the power of surface panel methods have established itera- 
tive inverse techniques for designing wing section geometry 
corresponding to a prescribed distribution of pressure (refer- 
ences 11-14). However, each of the existing three-dimensional 
wing design methods suffers from two or more of the following 
shortcomings : 

i) the calculations will not converge, 

ii) the designed geometry is unrealistically wavy chordwise 
and/or spanwise, 

iii) the designed geometry does not generate the prescribed 
pressure distribution near the wing leading edge, 

iv) the computing cost of each iteration cycle is 
substantial. 
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The new design method presented in this report was developed 
by McDonnell Aircraft Company (MCAIR) under contract to NASA, 
Langley Research Center. The method is believed to be the first 
that can reliably and efficiently solve the wing design problem of 
Figure 2 . The success of the method is attributed to the develop- 
ment of a formulation that overcomes the design dilemmas depicted 
in Figure 3 . 

The basic mathematical approach is similar to that of the 
two-dimensional design method developed for multi-element airfoils 
(references 8 and 15). However, the present design method has a 
significant new cost-savings feature that is especially appropri- 
ate for wing-fuselage configurations. The "perturbation analysis 
method" of references 9 and 10 is employed in each design itera- 
tion cycle. As will be demonstrated in Section 2, the perturba- 
tion analysis method is an accurate and extremely efficient method 
for analyzing large changes to wing section geometry. The present 
design method is designated the "perturbation design method" 
because it is an inverse formulation for the perturbation analysis 
method. The mathematical formulation and example design solutions 
are presented in Section 3 . 

Use of trade names or names of manufacturers in this report 
does not constitute an official endorsement of such products or 
manufacturers, either expressed or implied, by the National 
Aeronautics and Space Administrations. 
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2 . 


PERTURBATION ANALYSIS METHOD 


The standard solution approach to prescribed pressure design 
problems is to divide each iteration cycle into an analysis, 
pressure calculation step and an inverse, geometry correction 
step. In the two-dimensional method of references 8 and 15, an 
entire panel method solution is calculated during each analysis 
step. Furthermore, a new geometry- velocity perturbation matrix is 
calculated for each inverse step. In spite of the fact that the 
number of computations in each iteration cycle is a cubic function 
of the number of panels, the total computing cost is relatively 
small. The reason is that typical two-dimensional problems 
require fewer than one hundred panels. 

However, the number of panels required for wing-fuselage 
configurations is an order of magnitude greater, and the preceding 
design procedure would be extremely expensive. During the initial 
phase of the present contract, it became clear that a practical 
wing-on- fuselage design procedure could not be established on the 
basis of existing panel method technology. Consequently, the 
perturbation analysis method was developed. The method is an 
extremely efficient tool for analyzing the pressure distribution 
corresponding to a series of arbitrary, small perturbations to a 
baseline wing-fuselage geometry (Figure 4). The following 
features of the perturbation analysis method make it especially 
practical for application to an iterative wing section design 
method . 
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( 1 ) 


The computational expense for analyzing each successive 
geometry perturbation is more than an order of magnitude 
less than that of a conventional panel method, 

(2) The pressure distribution prediction accuracy is competi- 
tive with conventional surface panel methods for very 
large perturbations to wing section geometry, 

(3) A pre-calculated matrix of partial derivatives for the 
paneled baseline configuration is available. Each 
element of the matrix is the rate of change of potential 
at a boundary condition control point with respect to a 
geometry parameter perturbation. For design applica- 
tions, the geometry-potential perturbation matrix can be 
efficiently converted to a geometry-pressure perturba- 
tion matrix. 

The purpose of this section is to familiarize the reader with 
the fundamentals and power of the perturbation analysis method (a 
more complete description is available in references 9 and 10). 
An inverse formulation for the perturbation analysis method is 
reserved for the next section. The inverse formulation is the 
foundation for the iterative, prescribed pressure design method 
for wings and wings on fuselages. 


5 



2.1 MATHEMATICAL APPROACH - The perturbation analysis approach 
requires an initial baseline calculation of a conventional panel 
method solution for an arbitrary baseline configuration. 
Subsequently a matrix consisting of the partial derivatives of 
velocity potential with respect to geometry coordinates is 
calculated. The baseline solution and derivative matrix are 
calculated one time only and then stored for repetitive use. For 
each geometry perturbation, the solution surface distribution of 
velocity potential is constructed by multiplying the derivative 
matrix by a new right-hand-side. This procedure bypasses the two 
computationally expensive steps of a conventional panel method: 
calculating the influence coefficients and solving a large system 
of linear algebraic equations. 

Although the perturbation analysis method is appropriate for 
predicting the effect of arbitrary small changes to wing planform 
and fuselage geometry, the real power of the method is the 
accuracy with which large perturbations to wing thickness, camber, 
and twist can be analyzed. A simple two-dimensional example is 
the perturbation from a NACA 0012 airfoil to a circular cylinder 
illustrated in Figure 5. Unlike conventional "small disturbance" 
or "linearized" flow methods, predictions by the perturbation 
analysis method are accurate even at the leading edge of a wing . 
The reason is that only the surface potential - not velocity and 
pressure - is linearized with respect to geometry coordinate 
perturbations. It can be shown that the nonlinear terms are much 
smaller for potential than for either velocity or pressure. For 
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example, consider the flow about an elliptic cylinder. As 
illustrated in Figures 6 and 7, the variation of potential with 
thickness is exactly linear. 

Three computer programs are required in order to apply the 
perturbation analysis method (see Figure 8). The first program is 
the conventional MCAIR panel method for analyzing the baseline 
configuration. The second program employs an output solution file 
from the first program plus a differential mathematical formula- 
tion to calculate the matrix of partial derivatives of perturba- 
tion potential. For each baseline configuration, the first two 
programs generate an input file for the third program - the 
Perturbation Analysis Program. 

After the input file is generated, the first two programs are 
no longer required. The Perturbation Analysis Program can be 
executed repeatedly at low computing cost for the analyses of a 
series of perturbations to the panel corner coordinates (x,y,z)j 
of the complete aircraft configuration. The method used to 
analyze each perturbation is the same as the conventional panel 
method calculation with two significant exceptions. First, no 
influence coefficients are calculated; second, no large system of 
linear algebraic boundary condition equations is solved. Instead, 
the perturbation potential at each control point is calculated by 
linear extrapolation. The conversion of potential to surface 
velocity is based on numerical differentiation. 
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Usually the second program is executed twice for each base- 
line configuration, once at 0° incidence and once at 90° 
incidence. By employing the principle of linear superposition, 
the perturbation analysis is automatically performed at any inter- 
mediate angle of attack (Figure 9). 

2 . 2 EXAMPLE PERTURBATION ANALYSIS S OLUTIONS - Demonstration of 
prediction accuracy for the perturbation analysis method is pre- 
sented in Figures 10-12 . The baseline geometry is an isolated low 
aspect ratio wing with no camber and no twist. Coordinate 
displacements normal to the chord plane were applied in order to 
generate the fighter wing geometry illustrated in Figure 10b. The 
incompressible, potential flow solution for the fighter wing was 
calculated using both the conventional surface panel methods and 
the perturbation analysis method. The distributions of pressure, 
forces, and moment (Figures 12a and 12b) demonstrate that the 
perturbation analysis method is extremely accurate for very large 
perturbations . 

In order to evaluate prediction accuracy for a configuration 
with strong wing-fuselage interference, the paneled YAV-8B Harrier 
of Figure 13 was selected as the baseline configuration. The good 
accuracy of the conventional MCAIR panel method is demonstrated by 
comparing the calculated and experimental pressure distributions 
on the baseline wing (Figure 14). Compressibility effects were 
simulated by a Gothert correction, which is available in a later 
version of the panel method than the incompressible version of 
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reference 9. Approximately 420 seconds computing time on the CDC 
CYBER 176 were required for the analysis, where 537 panels on each 
side of the symmetry plane were used to model the YAV-8B. 

For application to the perturbation analysis method, the 
partial derivative of perturbation potential at each control 
point on the YAV-8B configuration was calculated with respect to 
the waterline displacement (zj) of each of the 150 panel corner 
points on the wing. The total computing time expended by the 
second program for 0° and 90° angles of attack was 1500 seconds on 
the CDC CYBER 176. 

The perturbation analysis method was then used to calculate 
the wing-fuselage pressure distribution corresponding to a large 
change in the wing geometry. The actual cambered supercritical 
wing with 8° twist (Figure 13) was perturbed to an uncambered, 
untwisted wing with constant NACA 0012 section geometry. It is 
obvious that this change is not recommended for the purpose of 
improving aerodynamic performance. It is simply intended to 
permit assessment of the accuracy of the perturbation analysis 
method for a relatively large change. As shown in Figure 15, the 
predicted pressure distribution agrees closely with the nearly 
exact potential flow solution calculated by the conventional 
surface panel method. However, the perturbation analysis required 
only eight seconds computing time compared to 420 seconds for the 
conventional analysis. On McDonnell Douglas computers, the total 
computing cost (including input/output) for the perturbation 
analysis was less than 1/20 the cost of the conventional analysis. 
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3. 


PERTURBATION DESIGN METHOD 


Consider a paneled representation of a wing or wing-on- 
fuselage at a fixed, arbitrary angle of attack. Suppose that 
prescribed pressure coefficients are assigned to the panel 
centers. The objective of the design method is to determine the 
change in wing section geometry - twist, camber, and thickness - 
that most nearly corresponds to the prescribed pressure 
distribution (Figure 2). 

In the preceding section, it was shown that the pressure 
distribution corresponding to large changes in wing section 
geometry can be accurately calculated by the perturbation analysis 
method. Conversely, it is reasoned that large changes in wing 
section geometry could be accurately designed by an inverse formu- 
lation for the perturbation analysis method. That approach is the 
basis for the perturbation design method of this report. The 
inverse formulation requires iteration because pressure is a non- 
linear function of the wing section geometry coordinates. 
However, the fact that perturbation potential is a nearly linear 
function eliminates the need to perform extensive computations in 
each iteration cycle. 

A schematic of the perturbation design method is presented in 
Figure 16 . The method can be used to solve very general aero- 
dynamic design problems. For example, the prescribed aerodynamic 
quantity at a panel center can be either a velocity component or 
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pressure coefficient. Arbitrary geometry parameters such as wing 
chord or fuselage shape can be selected for design. Furthermore, 
design constraints such as fixed camber or thickness can be 
imposed. Most design problems, however, are of the type 

illustrated in Figure 17. This type, designated the "standard 
wing design problem", is defined in detail below. The remainder 
of this section presents the mathematical formulation for the 
perturbation design method and example design solutions. 

3.1 STANDARD WING DESIGN PROBLEM - As illustrated in Figure 17, 
the region of panels subject to design is identified by corner 
points in the range ( iA' D a) £ (i#j) £ (iR/ js)/ where the limits 
( i&, ig, jj\, jg) are selected by the user. If and jg are points 
on the upper and lower surface trailing edge respectively, then 
the geometry of the complete wing section at each span station i 
will be designed. At the center of each panel in the design 
region, the desired pressure coefficient is prescribed by the 
user. 


The unknowns are Az at the panel corner points 
{(i-A'3A + D £ (i,j) £ (iB,jg-l)}. However, not all of the unknowns 
are permitted to be independent. As illustrated in Figure 18, 
less than one-half of the unknowns are independent. The remaining 
unknowns are generated by interpolation through the independent 
unknowns. On the span stations i = i^, ia+2, iA + 4# •••» ig* eac h 
dependent unknown Az(i,j) is established by least squares 
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quadratic interpolation through the path of points (j-3, 
j-1 , j, j+1, j+3) on the baseline configuration. For the 
remaining span stations (i = i^+1, i A +3, . .., ig-1), each unknown 
Az(i,j) is established by straight line interpolation through 
Az(i-l,j) and Az(i+l,j). For this type of interpolation to be 
accurate, it is implicitly assumed that the three points (i-l,j), 
(i,j), and (i+l,j) lie on nearly the same per cent chord line. 
Typical wing panelling is consistent with this assumption. A 
short computer program has been written to perform the 
interpolation on the baseline panelled configuration for all of 
the dependent variables. 

The reason for limiting the number of independent unknowns a 
priori is to prevent numerically unstable design calculations. 
Figures 19a and 19b illustrate two types of design instabilities 
that could occur if every panel corner in the design region were 
allowed to be an independent unknown. 

Consistent with the nomenclature of the perturbation analysis 
method (reference 9 ), each independent unknown perturbation is 
assigned an index kg. The value of Az for perturbation number kg 
is designated S^g. A schematic of the independent unknowns is 
presented in Figure 20. The objective of the perturbation design 
method is to calculate the values ( 1 _< kg _< NKS) that most 
nearly correspond to the prescribed pressure distribution. 
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3.2 MATHEMATICAL FORMUL ATION FOR PERTURBATION DESIGN METHOD - For 
any panelled baseline configuration, application of the perturba- 
tion design method requires that the arrays 4> i and 8<J> j_/9 Sj^g have 
been calculated a priori. <j> is the perturbation potential at the 
i^h control point and S^i/^S^g is the rate of change of 4>^ with 
respect to independent geometry perturbation number kg. The MCAIR 
3-D Subsonic Potential Flow Program and MCAIR 3-D Geometry 
Influence Coefficient Program will automatically calculate the 
required arrays and store them on a computer dish £ile. The 
perturbation design method can then be used to calculate the 
geometry perturbation that most nearly matches prescribed 
aerodynamic properties within the limitations of a minimal least 
square error. 


At any panel center selected by the user, one or more 
properties can be prescribed. The property can be either pressure 
coefficient (Cp) or velocity component in an arbitrary, specified 
direction. The prescribed value of an aerodynamic quantity at a 
panel center is designated Qp_^, where there is one index i for 
each prescribed value (1 < i < n DEs)* 


Arbitrary geometric constraints can be imposed upon the inde- 
pendent unknowns (S^ g ) . Each geometric constraint is expressed as 
a linear equation 



S W ] 


CHRS 


( 1 ) 
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Where n D , c ijy and CRHS are arbitrary values specified by the 
user. Depending upon the values specified, the constraint 
equation can be used to fix the cross-sectional area of a wing 
section, fix the thickness at one point, and so forth. 


The aerodynamic design problem can now be expressed in 

mathematical form. The objective is to calculate the array of 

independent geometric unknowns (1 _< kg <_ NKS) that will 

minimize the following function E. 

n DES / ) 


E - £ 1 WDES? • AREA. • (Q. - Q p ) 2 

i-1 ( 1 


n C0N 


r n E> 

- 

2 ) 

+ < 
j=l 1 

2 

WCON. . 

i 3 

11 

W -h q 
1 

Ci D ' Sk S (i D> - CRHS . 

J 


where (1) each i (1 <_ i _< n DES ) corresponds to one prescribed aero- 
dynamic quantity at one panel center. 


(2) each j (1 <_ j < hcqn) corresponds to a geometric 
constraint equation. 


(3) the weights WDESi and WCONj are specified by the user 
(typically, WDESi = i- 00 and WCONj » 1.00), 

(4) AREAi is the area on the baseline configuration of the 
panel corresponding to prescribed quantity Qp^, and 
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(5) (Q^ - Qp^) is difference between the calculated and 
prescribed values of an aerodynamic property. 

Iteration is used to solve for the array of unknowns S^g 
corresponding to the minimum value of E. As depicted in 
Figure 16, each iteration cycle is divided into an inverse step in 
which the geometry perturbation is calculated and a direct step in 
which the perturbed or updated geometry is analyzed. 


In the inverse step of each iteration cycle, the matrix of 
derivatives is calculated. Then the change in induced by a 
small perturbation to the array S^g can be expressed as 


NKS 

®i ■ £ [S-]- 

k =1 L9S k J 


By incorporating equation (3) in equation (2) and minimizing E 
with respect to ds^g, a system of linear, algebraic equations is 
established. The solution by standard matrix algebra provides the 
values dSfcg. The updated geometry is then analyzed by the 
perturbation analysis method. 

3Qi 

One might expect that recalculation of the matrix L r — — ] 

3Sk s 

during each iteration cycle would requre substantial computing 
expense. However, the following approach has proved to be both 
very efficient and accurate. 
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Which can be either 


Consider an aerodynamic quantity Q-^, 
pressure coefficient or a velocity component at the center of some 


panel. Figure 21 

illustrates 

the 

panel 

of interest 

and 

neighboring 

control 

points . 

If 

Qi 

is incompressible pressure 

coefficient 

, then 








Q i = 

C P 

1 - 

(V 2 

X 

2 

+ V + 

y 

v 2 ) 

Z 

(4) 


aQi 

-2(V X , 

V , V 

y 

z > * 

a 

(V ,V ,V ) 
v x' y' z 

(5) 


9S k 

K s 

as k 

K s 


If Qi is velocity component in an arbitrary, fixed direction (with 

direction cosines cosx, cosy, cosz), then 

-*• 

= V • (cosx,cosy,cosz) (6) 


aQ i 

as. 


- (cosx, cosy , cosz) 


as, 


( V V y'V 


(7) 


Equation (5) and (7) indicate that whether is pressure coeffi- 

aQi 


cient or velocity component, the problem of calculating 


3S k « 


is 


essentially reduced to calculating the velocity derivatives 


9V 


x 


av^ 


avr 


3Sk c ' 3S* ' 3S k , 


-) . This is discussed below. 


At the center of the panel of interest (Figure 21) , the veloc- 
ity vector V is 


V 


Ve+Ve+Ve 
xx y y z z 


V 

00 


+ V<p 


( 8 ) 
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On the basis of the perturbation potential $ at approximately 20 

neighboring control points j, V <J) can be calculated by numerical 

differentiation. This is expressed mathematically as 

^20 

■V* = E *j • (a j e x + b j e y + C j e z > (9) 

j=l 

where the scalars a j , b j , and Cj are functions of the corner point 
coordinates of the neighboring panels. A simplified analogue of 
equation (9) for one-dimensional numerical differentiation is pre- 
sented in Figure 22. (The three-dimensional numerical differentia- 
tion that is actually employed is based upon the assumption that 
neighboring panels stretch or shrink approximately the same per- 
centage. A more accurate algorthim could be introduced at the 
expense of greater compexity, but the present approach has been 
demonstrated to be adequate . ) 


By substituting equation (9) into equation (8) and 
differentiating, the desired velocity derivatives can be expressed 
as ^20 


as 


a 4> .. 

j=l ( s 


3 <W V Z > = S JsiJ- • (a j e * + b j e y + c j e ^ 


( 10 ) 


+ <J> . . 


(a .e + b .e + c .e ) 


3 3S. 3 X 3 y 3 z 


The only terra in equation (10) that requires substantial expense 
. 9<H 

to compute is . Fortunately, it is also the only term that is 

dSk S 
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nearly independent of perturbations to wing thickness, 
twist. Therefore, the precalculated baseline matrix 


camber, 

34>i 


and 


[ 


3 s k< 


1 that 


is available on computer disk file can be used in equation (10) 

during every iteration cycle. Substitution of equation (10) into 

, . _ x 3Qi 

equation (5) or (7) yields the desired value, -r— — . 

dS kc 


A significant feature of the preceding approach is that the 

3Qi 

accuracy of the calculated matrix — ] is competitive with an 

3S *S 

exact first order expansion during each iteration cycle. However, 
much less computing effort is required. 


In fact, the number of computations required for a complete 
iteration cycle is relatively small. The reason is apparent upon 
consideration of each calculation step in Figure 16. For example, 
consider the system of linear, algebraic equation to be solved for 
the perturbations to the independent unknowns. Typically, fewer 
than one hundred unknowns are sufficient for wing design, compared 
to several hundred for a conventional panel method analysis of a 
wing-fuselage. Also, consider the last step of each iteration 
cycle - analyzing the updated geometry. The extremely efficient 
perturbation analysis method is used for that calculation. 

The perturbation design method has been automated and is 
operational on the McDonnell Douglas CYBER 176. The computer 
program is designated "MCAIR Perturbation Design Program (Version 
1).“ A demonstration of the accuracy, efficiency, and numerical 
stability of the method is presented in the next section. 
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3.3 EXAMPLE DESIGN SOLUTIONS - Two example solutions by the 
perturbation design method are presented below. The calculations 
were performed on the McDonnell Douglas CDC CYBER 176. 

In the first design problem, the analytical pressure distri- 
bution for an infinite circular cylinder was prescribed. The 
baseline geometry is a NACA 0012 airfoil represented by 26 panels. 
The calculations converged to the minimum root-mean- square error 
in four iterations. As shown in Figure 23, the solution geometry 
is nearly circular. The slight discrepancy that does exist is 
attributable to the sparse panelling aft of 50% chord. The stream- 
wise geometry interpolation techniques of Figure 18 would have 
been substantially more accurate if 30 or more panels were used to 
represent the geometry, instead of only 26. 

In the second example, the objective is to design the fighter 
wing section geometry of Figure 10b by starting from the baseline 
wing panelling of Figure 10a. The fighter wing pressure distribu- 
tion at 0° angle of attack was prescribed at the center of the 208 
panels. The converged solution after 3 iterations is presented in 
Figures 24a and 24b. Although the designed geometry is close to 
the target, it is again apparent that slightly greater chordwise 
panel density would have permitted a better match on the lower 
surface at the tip. 
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Additional design examples have been successfully calculated 
using the panelling of Figure 10a as a baseline. A typical solu- 


tion requires 3 to 5 iteration cycles. The computing 
each complete iteration cycle is approximately 5 seconds, 
to 53 seconds for one conventional panel method analysis. 


time for 
compared 
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4. CONCLUSIONS 


The perturbation analysis and perturbation design methods are 
similar to classical thin wing theory in the sense that small 
disturbance "linearized" assumptions are employed. This mathe- 
matical simplification generates extensive computational savings 
for aerodynamic problems involving successive iteration, such as 
design. 

On the other hand, the restrictions of classical thin wing 
theory have been elminated. Compared to an exact potential flow 
solution, the present approach is quite accurate for thick wings, 
large leading edge radius or camber, and high angle of attack. 
The success of the perturbation design method is attributed to the 
inclusion of all significant first-order geometry-pressure peturba- 
tion terms in each iteration cycle. This leads to rapid solution 
convergence, in spite of the fact that the entire distribution of 
surface potential is constructed by simple linear extrapolation. 
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YAV-8B Surface Panel Modeling Upper Surface isobar Pattern 

Mach = 0.5 6.4* Angle-of-Attack 


Figure 1. Prediction Accuracy of MCAIR Panel Method 
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Given: • Fuselage Geometry 

• Baseline Wing 

• Angle*of- Attack 

• Prescribed Pressure Distribution (Cpj) 


Calculate: Wing Section Geometry 


Baseline 


^ 1 . J. 



\ ^ D< 


Desired 




Typical Wing Section 


Unknown Geometric 
Displacements bz 


x Prescribed Pressure Locations 


Figure 2. Objective of the Design Method 
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Figure 3. Common Barriers to a Successful Wing Design Method 




Objective 

• Subsonic Inviscid Analysis of Multiple Geometry 
Perturbations at Small Additional Cost 

Approach 

• Precalculated Baseline Matrix of Potential 
Derivatives ’"ayj 'Tjzjj 
Linear Extrapolation 

d0i 


0i + A0j)=<^i 




Figure 4. Perturbation Analytit Method 
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t Geometry 
Perturbation 
I (NACA-0012 
1 to Circle) 



■ ■■ Exact potential O Perturbation analysis 

flow solution method 

-3 


-2 


Cp — 1 


0 


i 


> 

) 



Figure 5. Simple Demonstration of Perturbation 
Analysis Method 

2-D Incompressible Flow 
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v = Voo+v<2> 


v “ =i — ^ T* 

I I L 



Thickness - t 


Figure 6. Variation of Flow Properties with Thickness 

Elliptical Cylinder at 0* incidence 
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Thickness — t 


Variation of Flow Properties with Thickness 

Elliptical Cylinder at 90* Incidence 









MCAIR 3-D 
Subsonic Potential 
Flow Program 

MCAIR 3-D 
fieomotry 

Influence 

Coefficient 

Program 

MCAIR 3-D 
Perturbation 
Analysis 
Program 

input 

e Baseline 
Geometry 

(x. y. z)j 

e Baseline 
Geometry 
(x, y. z)j 

e Baseline 
Potential 
Distribution 0j 

• Baseline 
Properties 
(x. y. z)j.0j 

• Derivative Matrix 
00j/9xj, d0j/dyj. 
d0j/dzj) 

• Geometry 
Perturbation 
(Ax, Ay. Az)j 

Approach 

e Conventional 
Panel Method 

• First-Order 
Expansion to 
Panel Method 
Formulation 

a Linear 
Extrapolation 
A0i« 

Z ((30j/oxi) Ax; 
j 1 ' 

+ (d0j/dyp Ayj 

+ (d0j/3zj) AZj] 

Output 

e Baseline 
Potential 
Distribution 0| 

e Baseline 
Aerodynamic 
Properties 

• Derivative 
Matrix 
(d0j/dxj, 

• Aerodynamic 
Properties 
of Perturbed 
Geometry 


Figure 8. Tire Thres Computer Program* for the 
Perturbation Analysis Method 




Figure 0. Linear Superposition of 0* end 90* Solutions 
Perturbation Analysis Method 
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Wing Section - NACA-0012 
208 Panels per Semi-Span 


Planform 

Aspect Ratio 2.9 

Taper Ratio 0.27 

Leading Edge Sweep 45* 

Trailing Edge Sweep 11.3* 




Figure 10 a. Baaellne Wing Paneling for Perturbation 
Analysis Test Casas 




i 

Root Section 

Tip Section 

Baseline Wing 

— . 

. __ 
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C 

Perturbed Wing 



(Fiflhter) 




Figure 10 b. Geometry Perturbation Tast Cate 
Fixed Planform Geometry 








Figure 11b. Baseline Fore* and Moment Distribution 
at 5* Angle-of-Attack 
MCAIR Panel Method 













Figure 12b. Fore* and Momant Distribution for Parturbsd Wing 
at 5° Angla-of-Attack 













Figure 14. Pressure Distribution on YAV4B 

6.4* Angle-of-Attsck Mach 0.50 
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Pressure at Panel Center Does Not Control 
Streamwise Slope Oscillations 


Figure 19a. Wavy Wall Design Instability 



Average Plane of Panel Does Not Control 
Corner Point Oscillations 


Figure 19b. FourComer Design Instability 
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Figure 20. Ordering of Independent Unknowns 
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Baseline (NACA 0012) 



Prescribed Pressure 

Converged Solution 
(4 Iterations) 


a) Pressure Distribution 



Baseline (NACA 0012) 
Target Geometry 

Converged Solution 
(4 Iterations) 


b) Geometry 


Figure 23. Circular Cylinder Design Solution 
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Planform 



Root 



Tip 


Baseline Wing Section (NACA 0012) 

Target Geometry (Fighter) 

o Converged Solution 

(3 Iterations) 


Baseline 


a) Geometry 


Prescribed Pressure 

Converged Solution 
(3 Iterations) 




JL 1 1— 1 1 I 

o 0-2 0.4 ae oj i.o 

x/c 



0.4 0.6 o.e 

x/c 


1.0 


b) Pressure Distribution 


Figure 24a. Example Wing Design Solution 
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